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SURFACE  IMPEDANCE  MODIFICATION  OF  PLATES 
IN  A  WATER-FILLED  WAVEGUIDE 


nmoDucnoN 

The  guidance  and  control  of  acoustic  wares  in  a  fluid  aediua  aaj  be 
accomplished  by  means  of  passire  materials  that  are  able  to  transmit,  reflect, 
and  absorb  acoustic  radiation  in  rarying  degrees.  There  are  limitations  to  the 
best  of  designs  using  passire  materials.  Actire  materials  offer  the 
possibility  of  expanding  the  range  of  control.  They  appear  especially  to  hold 
promise  at  low  frequencies,  where  passire  materials  may  lead  to  unwieldy  sise. 

During  the  pasi^  decade,  the  interest,  research,  and  industrial 
applications  of  actire  control  hare  increased  considerably,  especially  in  the 
field  of  actire  noise  control.  It  is  impossible  within  the  confines  of  this 
report  to  gire  eren  a  modest  introduction  to  the  relerant  literature. 
Therefore,  the  reader  is  referred  to  the  lists  of  references  in  the  articles 
giren  below,  and  especially  to  the  reference  bibliography  by  Guicking  [1]  . 
Fundamental  concepts  of  actire  sound  control  are  discussed  by  Ffowcs-Williams 
[2]  and  Scheuren  [3,4] . 

A  good  general  introduction  in  textbook  style  to  the  subject  of  actire 
control  of  sound  is  the  monograph  by  Nelson  and  Elliott  [5] .  It  gires  an 
extensire  discussion  of  interference  of  plane-ware  sound  fields  in  ware 
guides,  with  special  emphasis  on  the  question  of  where  the  power  in  the  field 
is  created  and  absorbed  (Chapter  5) .  In  connection  with  the  treatment  of  an 
absorbing  termination  at  the  end  of  a  duct  by  means  of  a  secondary  source,  the 
authors  mention  the  work  of  Bobber  and  others  at  the  Naral  Research 
Laboratory,  USED.  The  fundamental  properties  of  such  a  system  were 
inrestigated  by  Bobber  [8,7]  and  Beatty  [8] .  These  studies  were  applied  to 
test  sonar  transducers  under  traTeling-ware  conditions  in  a  tube  filled  with  a 
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liquid  at  high  pressure  [9] .  For  the  continuation  of  this  type  of  systea  in 
air,  Nelson  and  Elliott  refer  to  the  aork  of  Guicking  et  al.  [10,11].  Adaptive 
control  of  such  a  system  is  discussed  by  Orduna-Bustaaiante  and  Nelson  [12] . 

In  all  these  investigations,  it  is  typical  that  the  probing  of  the 
acoustical  field  is  accomplished  by  sensors  that  are  placed  in  the  field  sell 
aaay  from  the  actuator  surface.  In  contrast  to  this  arrangement  it  is 
attempted  in  the  present  study  to  influence  the  reflection  of  waves  by  a  plate 
in  a  water-filled  waveguide  by  means  of  layers  of  active  material  attached  to 
the  plate.  Thus,  the  complete  arrangement  of  sensor (s)  and  actuator  is 
designed  to  have  a  thickness  that  is  (much)  smaller  than  the  lateral  extension 
of  the  plate,  in  accordance  with  the  popular  notion  of  a  'smart  skin*. 

In  a  previous  study  [13] ,  a  no-transmission  experiment  in  a  water-filled 
waveguide  was  reported.  The  wave,  passing  through  a  single  layer  of  active 
material  (actuator) ,  is  suppressed  by  regulating  the  voltage  of  the  actuator 
through  a  feedback  loop  that  reduces  the  voltage  output  of  a  sensor  placed 
behind  the  plate  to  zero.  The  feedback  loop  is  closed  by  a  computer,  which 
performs  its  task  by  means  of  an  algorithm  from  complex-root-finding  concepts. 

Suppressing  the  wave  reflected  by  a  plate  is  more  involved,  since  one 
needs  two  items  of  information,  in  order  to  separate  the  amplitude  of  the 
reflected  wave  from  that  of  the  incoming  wave  [14] .  These  could  be  derived 
from  two  pressure  transducers,  or  from  one  pressure  transducer  and  one 
velocity  transducer.  It  might  be  expected  that  one  would  need  one  more 
transducer  yet,  acting  as  an  actuator,  in  order  to  produce  the  proper  acoustic 
wave  to  interfere  with  the  incoming  wave.  In  the  analysis  presented  here  it  is 
shown  that  one  may  establish  a  no-reflection  condition  by  means  of  only  two 
active  layers  attached  to  the  plate.  The  voltage  of  one  of  these,  the 
actuator,  is  governed  by  a  feedback  loop  based  on  the  same  algorithm  as 
referred  to  above,  which  uses  the  voltage  signal  from  the  other  layer,  the 
sensor . 

Two  aspects  may  be  distinguished:  first,  the  use  of  a  feedback  method 
based  on  the  technique  of  complex  root-finding;  and,  second,  a  mathematical 
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analysis  that  shows  how  a  combination  of  one  sensor  and  one  actuator  leads  to 
an  operable  feedback  loop  by  means  of  this  algorithm. 

After  a  description  of  the  equipment,  an  analysis  is  presented  to  show 
the  intended  operation  of  the  double  transducer.  Next,  experiments  are 
discussed  that  were  designed  to  show  the  feasibility  of  the  concept. 

lAVBGUIDS  AND  TRANSDUCERS 

A  sketch  of  the  waveguide  used  in  this  experiment  is  shown  in  Fig.  1. 
This  is  an  NEL-USRD  (Naval  Research  Laboratory,  Underwater  Sound  Reference 
Detachment)  type  G19  calibrator  [15] .  A  plane  wave  is  created  in  the  water- 
filled  tube  by  a  coil-driven  piston  in  the  bottom.  Reflection  of  this  wave  by 
the  free  surface  leads  to  a  standing-wave  pattern.  In  its  original  form,  the 
tube  wall  was  made  of  aluminum.  In  most  of  the  following  experiments,  a 
polymethylmethacrylate  (indicated  bf  the  trade  name  Incite)  tube  was  used. 
This  (transparent)  waveguide  was  originally  designed  for  laser  Doppler 
experiments . 


Fig.  1  -  BzpwiMntal  MrrugMvait  la  Q18  calibrator. 
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The  wave  speed  in  the  mediuv  contained  in  the  calibrator  (water)  is  not 

equal  to  the  wave  speed  in  a  aediun  of  unlimited  extent  due  to  the  elasticitj 

of  the  tube  wall.  To  establish  the  character  of  the  field  in  the  calibrator, 

the  pressure  near  the  piston  s  was  measured  hj  an  LCS  (Atlantic  Research) 

c 

miniature  hydrophone,  and  the  velocity  v  was  determined  from  the  response  of 

C 

an  accelerometer  mounted  on  the  side  of  the  piston  opposite  to  the  water 
column.  From  these  measurements,  one  nay  compute  the  specific  impedance  s^ 
seen  by  the  piston  by  s  =  p_/v  . 

The  results  are  shown  in  Fig.  2,  for  a  water  column  of  20.3  cm  height  in 
the  calibrator.  The  magnitude  of  the  impedance  is  given  as  a  function  of  the 
frequency.  The  solid  line  connects  the  experimental  points  indicated  by 
squares . 


I.WI  »  I  I  ■  III!  — •  -  '  - 

0.1  1.0  10. 

FREQUENCY  (kHz) 


Fig.  2  -  Spaeifie  iapadaBca  a  at  platos  svrfaea  aa  a  funetioa  of  fraqaanep- 

c 

Solid  earrat  eoapatad  froa  ■aaaarad  praoawa  and  Taloelty  at  piatoa. 

Daahad  ewrai  aeeordiag  to  ataading  vaTa  ralaa  a  >  e  taa  k  h, 

e  *^0  o  o 

wkara  e  la  araluatad  for  a  watar-fillad  infialta  alaatie  epliadar. 

9 

The  dashed  line  represents  the  theoretical  impedance  p  c  tan  k  h  for  a 

'  o  o  o 

water  column  of  height  h,  with  a  pressure  release  surface  at  one  end;  where 
is  the  density  of  the  water  and  s  propagation  speed  c^  is 

computed  from  a  model  for  axially  symmetric  propagation  of  waves  in  an 
infinite  water-filled  circular  elastic  cylinder  [16] .  The  ratio  of  outer  to 
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inner  radius  of  the  tube  is  1.125,  the  inner  radius  is  10.2  on,  the  height  of 
the  water  coluan  is  20.3  ca,  the  density  of  the  wall  aaterial  (Incite)  is  1200 

3 

kg/a  ,  its  Young’s  aodulus  is  4.74  GPa,  and  its  Poisson’s  ratio  is  0.316  (froa 
Ref. 17). 

Because  the  cross-sectional  areas  of  piston  and  tube  are  not  equal,  the 
theoretical  iapedance  was  noraalised  such  that  it  coincided  with  the 
experiaentally  aeasured  iapedance  at  the  low-frequency  end  (see  Fig.2) .  The 
factor  used  was  0.3,  not  far  froa  the  ratio  0.4  of  the  areas  of  the  piston  and 
tube  cross  sections.  One  sees  that  the  first  antiresonance  at  550  Hs  and  the 
first  resonance  at  1000  Hs  are  well  approxiaated  by  the  aodel.  At  higher 
frequencies  the  correspondence  between  experiaent  and  theory  appears  to  be 
lost. 


To  analyse  this  effect,  the  ware  speed  was  coaputed  froa  the  experiaental 
frequency  of  the  resonances  and  antiresonances  at  those  locations  where  they 
appeared  well  resolred.  This  is  coapared  with  the  coaputed  ware  speed  as  a 
function  of  frequency  shown  in  Fig.  3.  As  expected,  the  points  froa  the  first 
antiresonance  and  froa  the  first  resonance  are  close  to  the  theoretical  curre, 
but  the  points  froa  the  next  three  features  do  not  show  close  correspondence. 
The  experiaental  points  do  indicate  soae  decrease  of  ware  speed  with 
frequency,  but  they  do  not  confira  the  strong  dispersion  predicted  by  the 
aodel  above  about  1000  Hs.  The  cause  of  this  discrepancy  is  not  known. 
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fig.  S  o  Solid  ewToi  eoapotod  opood  o  of  oziollp  oyaMtrie  oovoo  io  o 

wotor-fillod  iafiaito  wovo  guido  wSth  olootie  epliadrieol  woll. 
Poiatoi  []  -  eoapotod  fro*  Mtiroooaoaoo  froqooaeioo  im  oacporiaostol 
eorro  (oolid  ewrro  ia  Pig.S). 

L  -  eoivotod  froa  roooaoaao  froqooMioo  ia  oaporiaaatel 
earro  (oolid  earvo  ia  Pig.S) . 


The  double  transducer  (Fig.  4)  is  constructed  fron  two  lasers  of  actire 
material,  each  3.3  an  thick.  The  active  material  is  NTK  Piesorubber  PR-306. 
('Piesorubber*  is  a  trademark  of  NTK  Technical  Division,  NGK  Spark  Plug  Co., 
Nagoja,  Japan.)  It  consists  of  PbTiO^  particles  embedded  in  a  neoprene 
elastomer  matrix.  The  center  electrode  is  common  to  both  transducer  disks,  and 
is  kept  at  ground  potential.  The  shields  of  the  transducers  and  the  shields  of 
the  coaxial  cables  are  electricallj  connected  together  and  to  ground.  The 
polarisation  of  the  two  transducers  is  antiparallel. 
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Fig.  4  -  Skateli  of  erooa  ooetioa  of  dooblo  troaodaeor. 

FEEDBACK  ALGOKITmi 

The  feedback  algorithm  used  in  this  experiment  is  more  extensively 
discussed  in  Ref. 13.  It  is  based  on  the  observation  that  the  purpose  of  the 
adaptive  feedback  arrangement  is  to  reduce  an  observable  quantity,  or  an 
expression  in  terms  of  observable  quantities,  to  zero.  In  the  case  of  a  no¬ 
transmission  condition  the  observable  quantity  is  simply  the  output  voltage  of 
a  sensor  placed  behind  the  plate,  the  transmission  of  which  should  be 
suppressed.  For  a  no-reflection  condition,  an  expression  in  terms  of 
observable  quantities  is  involved,  as  discussed  below.  In  any  case,  this 
expression  in  terms  of  observed  quantities  may  be  considered  as  the  dependent 
variable  w  of  an  analytic  complex  function  w.  The  independent  variable  of  this 
complex  function  is  the  input  voltage  to  the  actuator  z.  The  feedback 
algorithm  follows  from  the  notion  that  the  desired  result  is  equivalent  to 
finding  the  root  of  the  complex  function  w  =  f(z).  The  root  finder  used  in  the 
present  investigation  is  the  complex  equivalent  of  the  secant  method  in  real 
analysis . 


pnm  s.  pobbuaat 


In  this  root-finding  aethod  one  starts  with  two  pairs  of  rariables, 
and  (82**2^*  arbitrary  roltage  inputs  8^^  and  82  to  the 

actuator,  corresponding  to  two  ralues  w^^  and  W2  computed  from  the  obserred 
sensor  voltages.  A  new  actuator  voltag.;  8^  is  computed  by  82  = 

process  is  repeated  by  observing  the  corresponding  sensor 
voltage  and  computing  the  output  function  w^.  The  next  iteration  is  started  by 
setting  *2'^*2’  *3*^*2*  procedure  is  continued  until  the 

'cost  function*  Iwl  is  below  a  certain  preset  low  value. 

MODEL  FOR  OFERAIION  OF  DOUBLE  TRANSDUCER 

Unlike  the  suppression  of  a  transmitted  wave  described  in  Ref.  13,  it  is 
necessary  for  the  intended  suppression  of  a  reflected  wave  to  rely  upon  a 
mathematical  model  and  the  knowledge  of  certain  material  parameters  of  the 
active  material.  The  model  for  one-dimensional  wave  propagation  in  the  wave 
guide  used  here  is  discussed  in  detail  in  Appendix  A.  In  the  following 
discussion  the  definitions  and  results  of  Appendix  A  will  be  used. 

It  is  assumed  that  the  two  elastic  layers  are  made  of  the  same  material 
and  have  equal  thicknesses.  This  is  not  an  essential  restriction,  but  it 
simplifies  the  analysis  and  offers  better  insight  into  the  physics  of  the 
problem.  Results  for  the  relations  between  variables  for  the  general  case  of 
unequal  layer  parameters  may  be  found  in  Ref.  18.  The  model  shows  that  the  set 
of  eight  unknowns,  namely  the  forces  per  unit  area  and  the  velocities  at  the 
interfaces  at  both  sides  of  each  of  the  two  active  layers,  may  be  related  to 
the  controlled  or  observable  quantities:  the  voltages  of  actuator  and  sensor, 
and  the  current  through  the  actuator. 

The  algebra  gives  the  result,  not  surprising,  that  the  current  through 
the  actuator  is  mostly  due  to  the  capacitance  of  the  actuator  disk,  and  is 
little  affected  by  the  details  of  the  acoustic  impedances.  This  puts  excessive 
demands  on  the  accuracy  of  the  current  meter.  Therefore,  a  different  approach 
is  followed. 
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One  observes  the  two  voltaiges  and  of  actuator  and  sensor,  while 
there  is  no  external  voltage  impressed  on  the  actuator  and  the  terminals  are 
open-circuited.  Under  these  conditions  the  actuator  current  density  is 
known,  since  it  is  zero,  and  one  can  find  the  forces  and  velocities  at  the  two 
exterior  stirfaces  of  the  double  transducer.  Two  of  these  are  given  here. 


^21  ■  *11 


(^10'*~^20^  ^10~^20^  kd)/cos  kd 


2z 


13 


and 


^10^^20^  ~  ^10~^20^‘^°*  kd/ (1-cos  kd) 


'*^21  ■  2z 


13 


(1) 


(2) 


where  =  pc.  /(i  tan  kd) ;  p  is  the  density  and  c  the  propagation  speed  in 
the  transducer  material;  k  =  w/c,  d  the  transducer  thickness;  z^^^  ~  ^33/ * 
and  h  is  the  piezoelectric  constant. 

The  impedance  z^  presented  to  the  wave  by  the  plate  is  given  by  the  ratio 
of  force  and  velocity  at  the  interface  with  the  plate;  thus,  =  “^21^^21’ 
Then  z^  may  be  expressed  in  terms  of  a  non-dimensional  impedance  z^'  defined 
by  z^  =  z^'z^,  where  z^  is  related  to  the  parameters  of  the  layers  by 
z^  =  i  /7c  tan(kd/2),  in  the  following  form 


(Vj^+Vg)  cos  kd  +  (Vj^-Vg)  (1+cos  kd) 
p  ~  (1-cos  kd)  -  (Vj-Vg)  cos  kd 


z  = 


(3) 


For  small  values  of  kd  (this  condition  is  satisfied  in  the  experiments 


described  below) ,  the  expression  for  z  simplifies  to 

^  -3Vi^V2 

z„  =  0.5  i  upd  V  V  • 

P  ’'l  ■  ''2 


(4) 
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It  would  appear  reasonable  to  assume  that  this  impedance  is  constant  for 
a  long  period  of  time  (long  as  compared  with  the  time  needed  to  establish  the 
desired  no-reflection  condition) ,  and  thus,  once  computed  it  remains  arailable 
during  further  measurements  and  computations.  One  maj  obserwe  in  Bq.  (3)  that 
the  only  material  constants  needed  in  the  computation  of  the  plate  impedance 
are  the  propaigation  speed  of  plane  wares  (needed  for  computing  the  wave  number 
k)  and  the  density  of  the  plate  material  (in  the  definition  of  s^) .  When  kd  is 
small  Eq.  (4)  is  obtained,  and  only  the  pieaorubber  density  is  required  in 
addition  to  the  layer  thickness,  both  easily  measured  quantities. 

For  the  next  step  one  impresses  a  voltage  on  the  actuator.  Returning  to 
the  set  of  equations  A1-A6,  the  measured  s^  now  provides  a  relation  between 
force  and  velocity  at  the  plate.  Thus,  the  number  of  unknowns  reduces  by  one, 
and  may  be  left  unobserved.  The  algebra  again  gives  the  various  unknown 
quantities  expressed  in  terms  of  V2,  and  s^.  The  results  are  given  in 
Appendix  A.  The  expressions  are  rather  complicated,  but  it  suffices  to  notice 
that  the  surface  impedance  of  the  sensor  s^  =  ^12^^12  expressed  in 

terms  of  the  impressed  actuator  voltage  and  observed  sensor  voltage,  in 
addition  to  the  previously  determined  plate  impedance  s^. 

The  feedback  arrangement  reduces  here  to  the  problem  of  finding  an 
actuator  voltage  such  that  the  difference  between  the  value  of  t  and  a 
desired  impedance  is  reduced  to  zero.  For  a  no-reflection  condition  the 
value  of  z.  should  be  p  c,  where  the  subscript  sero  refers  to  the  values  of 
density  and  propagation  speed  of  the  medium,  which  in  an  elastic  waveguide  is 
affected  by  the  elasticity  of  the  wall.  The  relevant  function  w  may  be  written 
in  a  more  practical  form  by 

*  =  f j2"  *i^l2  ■ 

The  feedback  loop  is  designed  to  find  an  actuator  voltage  such  that 
the  absolute  value  of  this  function  w  (the  "cost  function*  in  control  theory 
language)  is  decreased  as  much  as  possible,  ideally  to  sero. 

In  principle  one  may  dial  any  desired  impedance  s^.  For  instance,  =  0 

would  impose  a  pressure  release  condition  and  s.  ■»  *,  equivalent  to  v.„  =  0, 

i  X2 

amounts  to  a  velocity  release  surface. 
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EZFBRIlffiNTS 

The  experiments  were  performed  according  to  the  analysis  presented 
before,  in  the  arrangement  sketched  in  Fig.l.  The  experiments  were  done  with 
and  without  a  backing  plate  behind  the  double  transducer.  The  operation  of  the 
feedback  loop  worked  satisfactorily.  Conrergence  to  the  desired  condition  was 
reached  in  about  three  steps.  Unfortunately,  an  important  criterion  for  the 
correctness  of  the  experiment  failed  to  be  satisfied. 

A  no-reflection  condition  established  at  the  plate  by  the  double 

transducer  necessarily  affects  the  impedance  seen  by  the  piston.  Specifically, 

the  impedance  of  the  water  column  at  the  piston  should  be  equal  to  the 

specific  acoustic  impedance  of  the  medium.  Notice  that  c  is  not  the  sound 

o 

speed  of  the  free  medium,  but  is  influenced  by  the  elasticity  of  the  walls  of 
the  calibrator,  as  discussed  above.  Moreover,  a  correction  factor  is  needed  to 
account  for  the  fact  that  the  piston  area  is  smaller  than  the  cross-sectional 
area  of  the  tube.  It  is  assumed  that  the  evanescent  waves  connected  with 
satisfying  the  boundary  condition  at  the  piston  have  little  effect.  The 
specific  impedance  at  the  piston  was  measured  by  means  of  a  miniature 
hydrophone  located  close  to  the  center  of  the  piston,  and  an  accelerometer 
mounted  on  the  back  side  of  the  piston.  In  all  cases  there  was  a  discrepancy 
of  at  least  an  order  of  magnitude  between  the  impedance  meaisured  directly  in 
this  fashion  and  the  impedance  p  c  entered  into  the  computer  program 
governing  the  feedback  loop. 

To  investigate  the  source  of  this  discrepancy,  the  pressure  and  velocity 
at  the  surface  of  the  double  transducer  (without  backing  plate)  oriented 
towards  the  free  water  surface  were  measured  by  the  miniature  hydrophone  and 
an  accelerometer  (Fig.  5} .  This  was  compared  with  the  values  computed  from  the 
observed  and  ¥2  of  the  actuator  and  sensor  when  no  external  voltage  was 
impressed  on  the  actuator.  A  typical  example  is  shown  in  Figs.  6  and  7.  It  was 
somewhat  encouraging  that  the  general  trend  in  the  directly  measured  values  of 
pressure  and  velocity  was  followed  by  the  values  computed  from  the  double 
transducer.  The  difference  in  the  magnitude  is  very  large,  though. 
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SlMteh  «f  satap  far  eo^ariaoa  «f  praaanra  aad  Taleeitj  eoaputad 
frea  Toltagar  of  dooblo  troaoduoar  (ao  earraat)  wltb  praoaara  frea 
alaiatara  bpdrepkoaa  aad  Taloeitp  froa  aeealareaabar. 


FREQUENCY  (kHz) 


Fig.  6  -  Coapariaoa  of  tha  praor’ira  ao  a  faaetioa  of  fraqnaaej  at  tha 

traaodaear  faea  oppooita  to  tha  iaeoaiag  wara  eoaputad  froa  Toltagaa 
ia  donbla  traaodaear  (aolid  liaa)  with  tha  aaaaaraaaat  by  a 
hydrophoaa  (daohad  liaa) . 
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fig.  7  -  Coapwrlaaa  of  tho  Toleeltf  m  a  foaetien  of  fTMiBaaej  at  tka 
traaodueor  faeo  oppooito  to  tho  iaeoaiag  waTo  eo^utad  from 
Toltagoo  ia  doablo  traaodaeor  (solid  lias)  witk  tka  ■aasaraaaat 
bp  aa  aeealaroaatar  (dasbad  liaa) . 

Extensive  study  was  Bade  of  the  possible  sources  for  this  discrepancy.  To 
exclude  possible  effects  of  propagation  in  the  walls  of  the  tube,  an  experiaent 
was  set  up  whereby  a  double  transducer  without  encapsulant  was  directly  placed 
on  a  shaker  in  air.  Even  in  this  simplified  case  the  ratio  of  values  for  the 
two  voltages  and  V2  could  not  be  reconciled  with  the  value  computed  from  Eq. 
(B4). 


It  was  concluded  at  last  that,  most  probably,  the  discrepancy  is  not  due 
to  any  errors  in  the  experiment  or  analysis.  Instead,  it  appears  that  the 
problem  is  due  to  deviations  in  the  actual  wave  action  in  the  transducer  from 
the  simple  one-dimensional  description  represented  in  the  model.  The  finite 
size  of  the  transducer  disks  allows  other  modes  than  just  a  pure  thickness 
vibration  to  enter  the  picture,  and  this  may  couple  into  the  voltage  measured 
across  the  transducer. 

In  Ref.  20  an  expe..'xz^w.t  with  a  single  piezoelectric  transducer  mounted  on 
a  shadcer  is  described.  The  authors  discuss  problems  with  discrepancies  between 
experiment  and  analysis  that  ao-e  due  to  deviations  from  one-dimensional  wave 
propagation.  They  found  that  very  careful  design  of  the  transducer  and  choice 
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of  a  specific  piesorubber  axe  required  to  obtain  results  that  concur  'with  the 
assuffled  one-dinensional  theory.  Further  discussions  with  one  of  the  authors 
(M.D.  McCollum)  supported  the  diagnosis  that  the  discrepancies  in  the  study 
reported  here  also  stem  from  deviations  from  one-dimensional  wave  theory. 

PI6UBE-QP-11ERIT  FOR  NO-BEFLECnON  EXPERIMENT 

In  Appendix  B,  a  derivation  is  given  for  a  figure-of -merit  (FOM)  in  the 
no-reflection  experiment.  It  is  assumed  that  the  transducer  layers  are  thin, 
i.e.,  kd  «  1.  For  both  steps  in  the  experiment  the  proper  expression  is  FOM  = 
d  hg^  /  ^33^*  Table  1  values  of  piesoelectric  coefficients  are  given, 
obtained  from  the  literature,  with  the  values  of  the  FOM  computed  for  d  =  3.175 
mm  (1/8*)  for  various  active  materials. 


Table  1.  FOM  for  no-reflection  experiment;  FOM  =  d  h^^  /  ^33^' 
assumed  that  the  thickness  d  is  3.175  mm  (1/8*). 


Material 

hgg  (GV/m) 

<=33° 

POM  (mV/Pa) 

PZT5 

2.15 

147. 

0.046 

BaTiOg 

1.56 

171. 

0.029 

PVDF 

0.97 

9.3 

0.33 

PR  305 

0.46 

2.4 

0.608 

PR  306 

0.38 

6.4 

0.188 

CONaUSIONS 

For  the  case  of  a  no-reflection  experiment,  the  feedback  algorithm,  based 
on  a  complex  root-finder,  works  properly. 

The  analysis  of  the  arrangement  of  two  active  layers  attached  to  a 
backing  plate  shows  that  it  is  theoretically  possible  to  establish  a  no¬ 
reflection  condition  by  means  of  feedback. 


14 


m.  nmoiAraoM  upoit  rsis 


The  attempts  to  experimentally  demonstrate  the  feasibility  of  this 
concept  were  inconclusive,  presumably  due  to  deviations  from  the  one¬ 
dimensional  wave-propagation  theory  implicit  in  the  model. 
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Appendix  A 

MODEL  FOR  NO-REFLECTION  EXPERIMENT 

BASIC  BQUAnOMS  AND  DSFINmONS 

In  Fig.  A1  the  geometry  of  a  plate  with  two  active  layers  is  shown. 


SOSOR  ACTlMItM  PUTE 


Fig.Al  Sksteh  of  intoroctlon  of  iBcomias  w»to  with  plots  oiwl  doablo 

tronoduesr.  Sorfoes  foreoo  p«r  nait  aroo  oro  iadieotsd  bj  > 
■urfocs  Tsloeitieo  bp  .  Th*  f Irot  oubaeripi  refsra  to  ths  aid*  of 

a  glToa  lapari  1  for  tha  aids  orimtad  towards  tha  < ng  wara,  3 
for  tha  oppoalta  sida.  Tha  aaeoad  sabaeript  rafara  to  tha  actira 
lapari  1  for  tha  actuator  and  3  for  tha  aanaor. 


A  harmonic  wave  with  presstire  p^  and  angular  frequency  v  impinges  on  the  sensor 
face.  The  pressure  of  the  reflected  wave  is  indicated  by  p^.  Assuming  thin- 
plate  theory,  for  which  the  thickness  of  the  plate  is  small  compared  with  the 
lateral  dimensions,  the  forces  per  unit  area  f^j  and  the  velocities  v^.  (i,j  = 
1,2)  are  related  by  the  following  set  of  equations  [19].  The  first  subscript 
refers  to  the  sides  of  a  given  layer:  1  for  the  side  oriented  towards  the 
incoming  wave,  2  for  the  opposite  side.  The  second  subscript  refers  to  the 
active  layers:  1  for  the  actuator  and  2  for  the  sensor. 
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^11  ■  *11  ’^ll  *  *12  ^21  *  *13  "^l 

^21  *  *12  ’ll  *  *11  ’21  *13  •^l 

^1  “  *13  ’ll  *  *13  ’21  *  *33  *^1 

^12  “  *11  ’l2  *12  ’22 

^22  *12  ’l2  *11  ’22 

^2  “  *13  ’l2  *  *13  ’22 

where 

*11  ~  f  (  ^  ^ 

*12  “  ^  ) 

*13  *  ^3  /  ^ 

*33  ~  /  (  i  **  ^33  ) 

p  -  density  of  solid 

C  =  (  Cgj”  / 

d  -  transducer  thickness 
k  =  w  /  c 

-  current  density  in  actuator 

-  voltage  across  actuator, 

^2  ~  voltage  across  sensor 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


with  the  following  definitions  of 
h  =  -(8T/8D)g  =  -(8B/8S)p 
e®  =  (8D/8B)g 

c®  =  (8T/8S)j, 


the  material  coefficients: 

D  -  electric  displacement 
B  -  electric  field 
S  -  strain 
T  -  stress. 
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The  boundarj  conditions  are 

£ 


12 
^12 


follows: 
P*  +  P  1 

p.  -  p 


^22  ’ 


and 


’ll  "  ’22 


(A7) 


where  s^  =  (p  ®  ^o®o*  *^o  influenced  bj  the  elasticity  of  the 

cylinder  wall. 

FIRST  STBP:  ND  VQLTACS  IMPIBSSED  ON  ACTUATOR 


In  this  case  the  current  density  is  sero.  One  solves  Eqs.  (A1-A6)  for 
the  variables  £221  fj^2»  ^22’  ’l2*  ^21’  ’21’  terms  of  the  two  voltages 
and  V2.  The  results  are 


-  <»11*»12)  O'l-V 

2*13 


f 


and 


'12 


2*l3('ir'X2) 


“  2.J3 


’21  " 


,  (»ll-«12)(Vi^V2)  -  »ll(^l-V 

2*i3(*ii"»12) 


(A8) 


To  simplify  these  results,  one  uses  the  following  properties  of  s^ ^ : 


*11  ”  *12  -  ^  (O.Skd)  , 

*11  *  *12  “  cot (O.Skd)  , 
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2  2  _  -  x2 

*11  ~  *12  ~  ’ 


*11 

*  *12 

1  *  COS  kd 

*11 

cos  kd  ’ 

*11 

-  cos  kd 

*11 

-  *12  ' 

1  -  CO*  kd  ' 

*11 

■"  *12 

-fl+cos  kd) 

*11 

II 

M 

1 

1  -  CO*  kd  * 

0-2  2  _  f  .2  2  co8(2kd) 

2*11  -  *12  -  1  -  cos(2kd) 


*ll[(^^V  ~  ^  kd)/cos  kd] 

2®.  A 


(Vl^V2)  ♦  (V1-V2)  CO*  kd/(l  -  CO*  kd) 

^3  ' 

*ll[^l‘^^2^  “*  kd)/co*  kd] 

2*«  n 


(AlO) 


(Vi^Y2)  -  (Vi-Vg)  cos  kd/(l  -  cos  kd) 

2*.  A 


The  surface  impedance  of  the  plate  s  is  found  by  z  =  -  fo-i/voi  and  the 

p  p  mX  mX 

input  impedance  to  the  sensor  surface  s^  by  s^  =  f^^  /  ^12' 


One  may  derire  a  compact  expression  for  the  non-dimensional  plate 
impedance  s^' ,  in  the  form 

*P  "  «ii  -  »i2  "  (Yi^V2)(1  -  CO*  kd)  -  (Yj-Vg)  cos  kd  ' 
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and  for  the  nondifflensional  sensor  impednnce  s  ’ 

®e  =  -  Zj2  "  (VJ+V2) (cos  kd  -  1)  -  (Vj-Vg)  cos  kd 


(A12) 


Notice  that  z  and  z  are  not  dependent  on  the  value  of  h»o,  but  only  on 

IP  6  ww 

the  values  of  p  and  c,  directly  in  the  factor  indirectly  through 

the  wave  number  kav/c. 


A  further  simplification  is  possible  when  kd  is  small,  namely 


.  3  V,  -  V2 

~  Y,  -  Y„  ' 


(A13) 


and 


•  -  3  V2 

*e  '  Y,  -  Y« 


(A14) 


The  dimensional  impedances  are  approximated  by 


Zp  -  0.5  i  up6. 


f-3Yj  -  Ygl 


-  ^2J 


(A15) 


and 


z^  -  0.5  i  w//d 


Ftr^J- 


(A16) 


Thus,  for  small  kd  the  results  for  z^  and  z^  depend  only  on  p  and  d, 
which  are  easily  determined  quantities. 


Notice  that  for  z^  =  0,  i.e.,  the  actuator  without  backing  plate  and 
exposed  to  air,  one  has  Y2  =  3Y^,  and  thus  z^  =  i  w/;(2d),  expressing  the  fact 
that,  for  small  kd,  the  input  impedance  is  just  the  inertia  of  the  double 
transducer  in  this  cue. 
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SBCOMD  STEP:  SQLUnON  WHEN  ■  IS  GIVEN 

P 

In  the  second  step  of  the  no-reflection  experiment  the  plate  impedance 
is  known  from  the  results  of  the  first  step.  One  enters  the  relations 

^11  "  ^22»  ^11  “  ■  ^22  *21  ■  ■  “p  ^21 

solves  for  the  variables  fj^gt  ^22’  ^22’  ^12*  ^21’  terms  of  Vj^,  Vg,  and 


2  2  2  2 

One  introduces  two  quantities  and  k^^  by  k^^  =  z^^  /*ii‘33 

k„^  -  z,,^/z-o*ao-  For  small  kd  one  has  k-^  -  k«^  -  k.^,  where  k.  is  the 
2  13  '  12  33  ^ 

thickness  coupling  constant,  k^  =  h^^  • 


The  results  are  expressed  in  compact  form,  by  again  using  z^'  ,  the  non- 


dimensional  impedBuice.  The  denominator  D  is  the  same  for  all  variables; 


=  fi-A.-y . [3--^ 


*13  ^*11'®12^ 


P  k^  1.2  ,2  ,2 


(A17) 


Only  the  expressions  for  f  ^12*  compute  the  input  impedance  to 

the  sensor  z^  =  f  12^^12’  given  here; 


*13^*ir*12^ 


,  f  2(l-l/k/)  (1  +  cos  kd) 

2  =  ""l^l  ^  V  "  "^2^  -  - 1  -  c^  kd - 


(l-l/kj  )  cos  kd  j  j 


p  1  -  cos  kd 


u  2  ,  2 


Dv-„  I  (1-1/k-^)  cos  kd 

^  -  .  -  CO.  M- 


(A18) 


f  1  11  '  fi  _2_1 

^"k2*u2  **p^'j^2 

r  1  ^^2  J  *1 J. 
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For  k(l«l,  and  -  k^  «  I,  one  has 


2  . 


Df-p  I  8+2z 


®13^*11‘*12^2  ^  P 


^  (kd  k^f  ’ 


(A19) 


and 


^12  * 

^12 -  ^  Y,(l+sJ  +  2Y. 


mW  -  *1 


•i3(*ir‘i2>  ■  ‘  ^ 


These  expressions  show  that,  although  Yl»Yn,  Y„  is  divided  by  the  small 
2  2  *  • 

quantities  (kd)  and  k^  ,  and  thus  the  tern  with  Y2  nay  be  of  the  sane  order 

of  magnitude,  or  larger  than  the  tern  with  Y. .  Thus,  it  is  necessary  to  know 
2  2 

(kd)  and  k^  to  sufficient  accuracy,  including  the  imaginary  part  of  k^. 
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Appendix  B 

FIGURE(S)-OF-MERIT  IN  NO-REFLECTION  EXPERIMENT 


In  order  to  derive  f igure(s)-of-Berit  (FOU)  one  requires  a  different 
solution  of  the  same  basic  equations,  where  now  the  voltages  and  Y2  are 
iinknovn.  These  solutions  may  also  serve  to  estimate  the  expected  or  necessary 
voltages  for  sensor  and  actuator,  respectively,  to  predict  the  measurement 
values  in  a  given  experiment. 


In  the  first  part  of  the  experiment  one  determines  from  the  measured 
Yj^  and  Yg,  while  =  0.  One  solves  Eqs.  (A1-A6)  for  Yj^,  Yg,  ^22'  ^22’  ^21' 
and  V221  terms  of  f^2i  Introducing  a  non-dimensional  impedance  z^  by 


I  _ 


=  »«  /  (*ii  ■  *19)- 


12-* 


and 


^12*13 


’ll 


1  +  z 


2(1  +  Zp')  cos  kd  +  (2  -  Zp'/cos  kd) 


(Bl) 


^2  = 


^12*13 


’ll 


2  cos  kd  (1  +  z  ')  +  (1  -  z  ') 

_ ^ _ EJ _  P 

2(1  +  Zp')  cos  kd  +  (2  -  Zp'/cos  kd) 


(B2) 


To  reduce  the  signal-to-noise  ratio  the  voltages  should  be  as  large  as 
possible.  One  sees  that  for  a  given  incoming  wave,  represented  by  fj^2>  '^Be  two 
voltages  will  be  proportional  to  the  ratio  *Bich  therefore  may  serve 
as  a  FOM. 


For  small  values  of  kd  this  FOM  simplifies  to 


FOM 


^33! 

D  • 
®33 


(B3) 


Thus,  for  small  kd  one  may  write: 
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1  +  * 


=  FOM 


4  s 


3  *  g 


s  FOM 


4  +  «. 


The  ratio  V„/V-  Taries  from  3  for  s  '=0  to  1  for 

2  l  P  P 

Not  surprisingly  (see  Ref .  22) ,  the  other  variables  may  he  expressed  as 
f  olloirs : 

^11  ^22  i  PC  tan  kd  +  s 

^  ,>c  »  i  t»n  ka  • 


For  small  kd,  this  reduces  to  the  impedance  i/>tfd  (  1  s^'  /2) 


Similarly, 


i  pc  tan  2kd  s 


*e  ”  Vj2  ~  pc  *  iSp  tan  2kd 


Further, 


1  +  cos  kd(l+a  ) 


^12  2  cos  kd  (!■*■  cos  kd)  ^  s  (2  cos^kd-1) 


For  small  kd  this  reduces  to 


*22  *  * 


'12  A  *  g 


The  ratios  of  the  other  velocities  to  V22  8^8 

Y  I 

=  [cos  kd  +  s  (cos  kd-l)]~^, 
»oo  P 
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I  2 

t^2  kd  (1-cos  kd)  1  -  2cos  kd 

_  _  j  ^ 

22  cos  kd  *  8^ (cos  kd-1) 

For  small  kd  these  ratios  reduce  to  1  and  -1,  respectively. 


(BIO) 


In  the  second  part  of  the  no-reflection  experiment  the  value  of  z  is 

P 

given,  and  the  value  of  the  voltage  on  the  actuator  is  adjusted  to  reach  the 
desired  z  .  Again,  in  order  to  improve  the  signal-to-noise  ratio  one  wants  to 
optimize  the  voltage  V2,  given  ^Q'^roducing  a  non-dimensional  input 

impedance  by  zj  -  “  *12^'  finds  by  solving  the  Eqs.(Al-A6)  for 

^2’  ^22’  ^22*  ^21’  *^1’  “  of  fj* 


^2  = 


*13  ^12  *e'  ■  ^ 


(Bll) 


Thus,  one  finds  here  for  a  FOU  the  ratio  z^^^  /  ^12*  *^foh,  for  small  kd, 
reduces  to  the  same  value  for  the  FOU  as  in  the  first  part  of  the  experiment. 
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